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Abstract - Using infrared spectroscopy, we investigate bottom gated ABA-stacked trilayer 
graphene subject to an additional environment-induced p-type doping. We find that the 
Slonczewski-Weiss-McClure tight-binding model and the Kubo formula reproduce the gate voltage- 
modulated reflectivity spectra very accurately. This allows us to determine the charge densities 
and the potentials of the 7r-band electrons on all graphene layers separately and to extract the 
interlayer permittivity due to higher energy bands. 



a 

O 
o 



> 

oo 

m 
O 

ON 

O 

(N 



Introduction. — Trilayer graphene has recently at- 
tracted much interest because its electronic structure is 
distinctly different from the bands found in more stud- 
,ed monolayer and bilayer graphene. Interestingly, the 
two major stacking types of trilayer graphene, namely the 
Bernal (or ABA-)stacking and the rhombohedral (ABC- 
)stacking, result in rather different band properties [1]. In 
rhombohedral graphene, a sizeable bandgap can be gener- 
ated and tuned with an external electric field [2-5], while 
in Bernal stacked graphene the dominant effect is gate- 
tunable band overlap [2,6-11], resulting in a semimetal- 
lic behaviour. Trilayer graphene presents the simplest 
case where the complete set of Slonczewski-Weiss-McClure 
(SWM) tight-binding parameters are needed to describe 
the band structure, as it is the case for three-dimensional 
graphite [12, 13]. However, in contrast to graphite, tri- 
layer graphene allows exploring different doping regimes 
and the effect of an electric field on the band structure. 
Thus it might be beneficial for understanding the complex 
graphite physics to study its trilayer counterpart. 
In this work we focus on the ABA-stacked trilayer 
graphene, the form most often found in nature, which is 
plane symmetric with respect to the middle layer. Our 
goal is to determine experimentally optical properties as a 
function of an external electric field and check the ability 
of the existing theories to reproduce them quantitatively. 
The sample is a bottom gated exfoliated flake. Its charge 
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Fig. 1: (a) Schematic representation of bottom gated three- 
layer graphene used in this study. The layer numbering used 
in the text is given, (b) The resistance of the sample as a 
function of the gate voltage at 20 K. 



neutrality point is strongly shifted with respect to zero 
gate voltage, which allows us to explore the high doping 
regime not achieved so far in a combination with infrared 
spectroscopy. By treating the electrostatic potentials of 
each of the three layers as gate voltage-dependent free pa- 
rameters, we achieve a very good match between the ex- 
periment and the tight-binding model. With the help of 
a basic electrostatic model for trilayer graphene [7, 8] we 
find an accurate value for the static perpendicular inter- 
plane permittivity e r due to the cumulative polarizability 
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Fig. 2: Reflectivity of trilayer graphene at 20 K for different 
gate voltages of -80, and +80 V. 

of all electronic shells except the low-energy 7r-bands, con- 
sidered explicitly by the tight-binding model. 

Experimental results. — A large, several tens of mi- 
crons, flake of trilayer graphene was obtained by mechani- 
cal exfoliation [14] and deposited on a Si substrate covered 
by a 300 nm layer of silicon oxide (Fig. 1(a)). Electrical 
contacts were deposited by evaporating gold through a 
shadow mask. The stacking was identified as of the ABA- 
type based on the optical spectra as it is discussed below. 
The resistivity curve (Fig. 1(b)) shows that the sample 
is strongly hole doped so that the charge neutrality point 
cannot be achieved even by applying V g = +80 V to the 
gate. Annealing the sample several times at 200 °C in a 
H2 - N2 atmosphere did not change significantly the resis- 
tivity curve. We therefore attribute the extrinsic doping 
either to gold atoms, generally known to dope graphene 
positively [15, 16], scattered below the shadow mask onto 
the graphene surface during the contact evaporation or to 
a layer of water or hydrocarbon contamination trapped 
below the flake [17]. Although not being able to reach the 
electron-type doping, we benefited from this situation by 
measuring in the regime of high hole doping (~ 3 • 10 13 
cm -2 ), which would not be achievable with an initially 
undoped sample. 

The reflectivity measurements at various gate voltages be- 
tween -80 V and +80 V with a step of 20 V were performed 
at 20 K in a helium-flow cryostat mounted on an infrared 
microscope coupled to a Fourier transform infrared spec- 
trometer. The infrared spot of about 15x15 /im was kept 
far from the flake edges and electrical contacts. Each mea- 
surement on graphene was followed by recording reference 
spectra on the bare substrate and on gold. A black-body 
Globar light source was used between 0.1 and 0.7 eV and 
a tungsten lamp was employed between 0.4 and 1.5 eV. 
The reflectivity of trilayer graphene normalized by the 
bare substrate is plotted in Fig. 2 for gate voltages of -80, 
and +80 V. Although the changes between the spectra 
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Fig. 3: Experimental differential reflectivity for pairs of con- 
secutive gate voltages. The value given for each curve is the 
average of the two voltages (separated by 20 V) used to gen- 
erate the difference. For clarity each curve is shifted and the 
zero levels are indicated by dotted lines. The points represent 
the experimental data while the solid red line is the result of 
the fit. 

are caused by the modification of the optical properties of 
graphene, some profound spectral features present in all 
spectra are due to the substrate. For example, the strong 
structure at about 0.1 eV is due to optical phonons in 
SiC>2 • The deep minimum at 0.7 eV, where the presence 
of graphene reduces the reflected signal by almost 40 per- 
cent, is due to the Fabry-Perot interference in the oxide 
layer. 

Even though the substrate somewhat complicates the in- 
terpretation of the measured spectra, some qualitative ob- 
servations can be made without data modeling. First, at 
energies below about 1 eV the reflectivity systematically 
decreases as the gate voltage increases. This is directly 
related to the decreasing metallicity and Drude weight as 
the concentration of holes is reduced. Second, above 1 eV 
the curves essentially overlay showing that the electric field 
has little or no influence on the optical properties at high 
energies. Third, in the region between about 0.4 and 0.8 
eV, a more complex gate voltage dependence is observed, 
which is due to interband transitions in graphene. 
In order to emphasize the effect of the gate and to extract 
more information from the data, we show in Fig. 3 the 
differential (or contrast) spectra for all pairs of consecutive 
gate voltages. They are defined as follows [18]: 

AR _ A V g + 10V) - R^j, V g - WV) 

R [ ' 9> R(u, V g + 10V) + R(u, V g - 107) ' 1 ' 
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The curves are marked by pronounced spectral structures 
showing a strong gate dependence. For example, the posi- 
tion of the dip at around 0.6 eV shifts by almost 0.1 eV as 
the gate voltage is swept from +70 V to -70 V. The shape 
of the structures is also strongly influenced by the gate. 
Such a non-trivial behavior of the contrast spectra is due 
to the fact that not only the chemical potential is modified 
by the gate voltage but also the band structure is affected 
by the gate-induced electric field. Thus quantitative mod- 
eling is required in order to identify spectral features and 
to get further insight into the effect of the gate on the 
electronic structure of the sample. Below we discuss this 
analysis in detail. Looking ahead, we note here that a very 
good agreement between the data and the tight-binding 
model of the ABA-stacked graphene is achieved (red solid 
lines in Fig. 3). Given that the band structures of ABC- 
and ABA-stacked graphene are entirely different [1], we 
conclude that the sample is ABA-stacked. 

Optical data modeling. — To model the optical data 
we use a tight-binding description of trilayer graphene [7,8] 
based on the SWM-model of graphite [12, 13]. The hop- 
ping terms 70, 71, 72, 73, 74 and 75 are illustrated in Fig. 
4. Additionally, A represents the difference between the 
on-site energy for the two sublattices A and B (marked 
with red and blue colors respectively). To account for the 
effect of an external electric field we add potentials U\,U2 
and U3 to the first, second and the third layers respectively 
[7]. In the basis of |^) , \B{) , \A 2 ) , \B 2 ) , \A 3 ) and \B 3 ) 
the Hamiltonian of ABA-stacked trilayer graphene reads: 
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where k is the momentum with respect to the K-point of 
the Brillouin zone and the vectors d 1^,3 connect the near- 
est neighbor atoms as shown in Fig. 4. 
In order to obtain the optical properties of this system, 
within the linear response theory, we use the Kubo for- 
mula for the optical sheet conductivity a(u). We calcu- 
late the Drude and the interband terms using the following 
expressions 
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Fig. 4: The crystal lattice of ABA-stacked trilayer graphene 
with the definition of the hopping parameters. 
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where a = e 2 /Ah w 6.08 • 10~ 5 f2 _1 is the universal AC 
conductance of monolayer graphene [19], is the band 

energy (i = 1,..,6 ), /(e) = (l + exp(f^)) is the 

Fermi-Dirac distribution and /x is the chemical potential. 
F and r^) are the electronic broadening parameters for 
the interband and the intraband transitions respectively. 
For simplicity, the interband broadening is taken to be 
momentum-, energy-, band- and doping-independent. 
The total charge density n to tai and the charge densities of 
each layer, m (Z = 1,2,3), can also be calculated through 
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where we subtract \ in order to count the doping level 
with respect to half filling of the n bands. ipA l and ipB t 
are the projections of the total electronic eigenfunction on 
the sub- lattices Ai and £>;. 

Once the conductivity is determined we can calculate the 
optical reflectivity of the complete graphene/Si02 /Si sys- 
tem using Fresnel eq.s and consequently compute the dif- 
ferential reflectivity (Equation 1). As AR/R values are 
rather small it is useful to introduce the so-called sensitiv- 
ity functions [18] /3i(w), /32(w), which are determined by 
the substrate, and to employ the linearized formula 
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When fitting the experimental data, we fixed the SWM 
parameters to the values published in Ref. [20], which we 
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Table 1: SWM band parameters used in this work to fit the 
optical spectra. The values are taken from Ref. [20]. 



Parameter 


Value (eV) 


7o 


3.12 


7i 


0.377 


72 


-0.0206 


73 


0.29 


74 


0.12 


75 


0.025 


A 


0.0138 



reproduce in Table 1 for the sake of convenience. We also 
fixed the Drude broadening to 5 meV [21] as its value 
has almost no influence on the optical data in the stud- 
ied spectral range. In contrast, the interband broaden- 
ing T strongly affects the data, resulting in the overall 
smoothening of the spectral structures. It was therefore 
taken as an adjustable parameter, for which the best value 
was found to be 45 meV. 

The most important detail of our fitting procedure is that 
we treat the chemical potential and the plane potentials 
Ui as independent adjustable parameters at every value 
of V g , in order to extract the actual gate voltage depen- 
dence of these quantities from the optical data. It is easy 
to see that shifting the values of fi, U\, U% and U3 by 
the same amount does not change any physical properties 
such as the optical spectra and the charge density. To re- 
move this ambiguity, we imposed an additional condition 
U 1 + U 2 + U 3 = 0. 

The best fits are shown in Fig. 3 with the solid red lines. 
One can notice an excellent match between the experiment 
and the model, including all spectral features at all gate 
voltages. It is worth to mention that although the differen- 
tial reflectivity between 0.15 and 0.4 eV is rather feature- 
less it is distinctly non-zero. This is caused by the gate 
variable Drude contribution and the model reproduces this 
fact very well. 

As all the curves were fitted simultaneously, the obtained 
parameters were well defined and converged to the same 
numbers, which we carefully tested by using different fit- 
ting procedures and by varying the initial parameter val- 
ues. The robustness of these results allows us to develop a 
deeper analysis of the response of graphene to the external 
electric field which we present in the next section. 

Discussion. — The optical conductivity spectra cor- 
responding to the fits described above are shown in Fig. 
5(a) for each value of the gate voltage. For illustration pur- 
poses, the curves are normalized by 3 ■ erg. At low photon 
energies but above the Drude peak, the optical conductiv- 
ity is significantly below 3 • <tq , while at high photon ener- 
gies it approaches this value. This is related to the Pauli 
blocking, which results, in the case of monolayer graphene, 
in a conductivity step at 2fj, [22]. Additionally, significant 
spectral features are the peaks at approximately 0.25 eV, 



0.42 eV and 0.6 eV, which are denoted respectively as A, B 
and C in Fig. 5(a). They can be assigned to transitions be- 
tween almost parallel bands, which are illustrated in Fig. 
6. The energy of the strongest peak (C) is close to V%Ji, 
in agreement with the theoretical predictions [6, 23, 24] 
and experimental data [4,11,25,26]. The peaks noticeably 
shift as a function of the gate voltage. For example, peak 
C moves from 0.55 eV at V g = +80 V to 0.6 eV at V g = 
-80 V. As can be seen in Fig. 6, this effect stems from 
the field-induced displacement of the bands with respect 
to each other. The transitions are broad because the scat- 
tering parameter T is large, as compared to the numbers 
obtained in low doped samples [4]. We believe that extra 
scattering present here is due to the charged impurities 
that also give rises to the elevated doping. 
Let us now inspect the dependence of the doping and the 
intcrplane potentials on the gate voltage. In Fig. 5(b) we 
plot the chemical potential as a function of the applied 
gate. It is always negative with a high absolute value and 
moves gradually towards the charge neutrality as the gate 
voltage is increased. This corroborates our initial obser- 
vation of a high level of hole doping. 

The total charge carrier density ntotai obtained using 
Equation 6 is shown in Fig. 5(c). At zero voltage the 
concentration of holes is about 3 • 10 13 cm -2 . It changes 
linearly with V g with a slope dn to tai/ dV 9 equal to 7.0 • 
10 10 cm~ 2 V _1 . Remarkably, this value, which is based 
on the optical data only, agrees extremely well with the 
expected slope of 7.2 • 10 10 cm~ 2 V _1 for the actual thick- 
ness of the oxide layer [14], thus providing extra credence 
to our analysis. From this dependence, one expects to 
reach the charge neutrality point at V g « +450 V, which 
is not achievable in the present experiment. 
In the same Fig., we present the charge density on each 
layer separately according to Equation 7. As expected, the 
layer closest to the gate shows the strongest charge vari- 
ation as a function of the external field, while the effect 
of the gate decreases strongly with increasing layer index 
due to internal screening. This is in good agreement with 
theoretical predictions [8] . 

In Fig. 5(b) we plot the potential difference between the 
outer layers, Ai = XJ\ — U3, and twice the difference be- 
tween the average of the potentials of outer layers and the 
one of the central layer, A 2 = U\ + U3 — 2C/ 2 , introduced 
in Ref. [7] . These parameters fully characterize the distri- 
bution of the potential in the sample. We first note that 
the value of Ai is not zero but 0.25 eV in the absence of 
the gate voltage. Such a large potential asymmetry can 
be attributed to environmental impurities in contact with 
the sample, whose charging effect is not symmetrical with 
respect to the graphene plane. Since Ai depends on V g , 
the gate also contributes to the total electric field across 
the sample. The physical meaning of the parameter A2, 
which is also finite at V g = and shows a noticeable gate 
voltage dependence, will become clear in the following. 
Within a simple electrostatic model [7, 8, 27] where each 
graphene layer is treated as an infinitely thin charged 
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Fig. 5: (a) The calculated optical sheet conductivity at different gate voltages obtained from the fits shown in Fig. 3 as described 
in the text. The arrows mark the dominant peaks corresponding to the interband transitions indicated in Fig. 6. (b) The 
interplane potentials Ai, A2 and the chemical potential /1 as a function of the gate voltage, directly obtained from the optical 
fits. Also shown are the derived values of the external potential Ai xt and the interspace permittivity e r . (c) The charge densities 
ri; of each layer and the total density n to tai as a function of the gate voltage. 



plane the electrostatic potentials can be straightforwardly 
related to the charge densities: 

A! = — (m-naJ + Ar*, (9) 

EoS r 

and 

e 2 d 

A 2 = -—n 2 . (10) 

where e is the elementary charge, £0 the vacuum per- 
mittivity, d w 0.335 nm the interlayer distance and e r 
the interplane permittivity between graphene layers due 
to the screening effect of all electrons except the 7r-band 
electrons. The quantity Af xt takes into account the ef- 
fect of the gate voltage and the environmental charges, 
which would produce a potential difference exactly equal 
to Af xt over the distance of 2d in the absence of electro- 
static screening by graphenc. 

Interestingly, all the quantities that enter Equation 10, ex- 
cept e r , are known from the infrared data, which allows us 
to derive the interplane permittivity. From Fig. 5(b) we 
cannot identify any dependence of e r on the gate voltage, 
and its value is 3.7 ± 0.1. 

To our knowledge, there are no direct measurements of 
this quantity, which is however an important parameter 
needed to calculate the electronic screening in graphene. 
One should note that no relation analogous to Equation 
10 exists in the case of bilayer graphene allowing such a 



straightforward determination of e r . It turns out that the 
high doping level of our sample is beneficial for the accu- 
racy of this method since e r is proportional to the ratio of 
712 and A2 , which would both be very small in an undoped 
sample. 

Once the value of e r is determined, we can calculate 
Af xt with Equation 9 for each gate voltage (light gray 
squares in Fig. 5(b)). Af xt decreases linearly as a func- 
tion of V g with a slope of about -0.0014 eV/V. Within 
the same electrostatic model, the slope is expected to be 
dAfVdVg = ~{dee r )-(e b /L b ) = -0.0012 eV/V, where e b 
= 3.9 and L b = 300 nm are the dielectric constant and the 
thickness of the oxide layer respectively. The experimental 
and the theoretical values are thus very close, supporting 
the model. 

It is interesting to follow the modification of the band 
structure by the gate voltage. The bands corresponding 
to the best fits of the optical spectra are displayed for V g 
= -80 V and V g = +80 V in Fig. 6 by solid lines. Al- 
though the bands for the presented voltages are qualita- 
tively similar, their positions and shapes, especially near 
the K point, are significantly influenced by the electric 
field. In both cases, there is a gap of about 50 meV be- 
tween the valence and the conduction bands. However, 
the chemical potential does not lie inside the gap, giving 
rise to a metallic, rather than semiconducting electronic 
properties. Notably in our highly doped sample the gap is 
much larger than the prediction for weakly doped ABA- 
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Fig. 6: Electronic bands at V g = -80 V and + 80V calculated 
based on the fits in Fig. 3. The horizontal lines indicate the 
chemical potential. The arrows correspond to the optical peaks 
indicated in Fig. 5(a). 



stacked trilayers [8,28,29]. 

In conclusion, we investigated highly p-type doped ABA- 
stacked trilayer graphene with infrared spectroscopy. By 
fitting the experimental reflectivity spectra with a tight- 
binding model and the Kubo formula we extract the in- 
trinsic parameters such as the carrier density and poten- 
tials of individual graphene layers. With a simple elec- 
trostatic model we can account for the effect of environ- 
mental charging and determine the interlayer permittiv- 
ity between two graphene layers due to all the electronic 
shells except the 7r-bands. We show that, due to a large 
potential difference across the layers, a sizeable band gap 
in ABA-trilaycr graphene is induced. Our work demon- 
strates new abilities of infrared spectroscopy in studying 
electronic screening in graphene multilayers. 
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